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ABSTRACT 

The two-point clustering of dark matter halos is influenced by halo properties besides 
mass, a phenomenon referred to as halo assembly bias. Using the depth of the gravi¬ 
tational potential well, Umax, as our secondary halo property, in this paper we present 
the first study of the scale-dependence assembly bias. In the large-scale linear regime, 
r > 10Mpc//i, our findings are in keeping with previous results. In particular, at the 
low-mass end (Mvir < M^ow « halos with high-14nax show stronger 

large-scale clustering relative to halos with low-Umax of the same mass; this trend weak¬ 
ens and reverses for M^ir ^ Mcoii- In the nonlinear regime, assembly bias in low-mass 
halos exhibits a pronounced scale-dependent “bump” at 500kpc//i — 5Mpc//i, a new 
result. This feature weakens and eventually vanishes for halos of higher mass. We show 
that this scale-dependent signature can primarily be attributed to a special subpopula¬ 
tion of ejected halos, defined as present-day host halos that were previously members 
of a higher-mass halo at some point in their past history. A corollary of our results is 
that galaxy clustering on scales of r ~ 1 — 2Mpc//i can be impacted by up to ~ 15% 
by the choice of the halo property used in the halo model, even for stellar mass-limited 
samples. 


1 INTRODUCTION 


The halo model provides a connection between dark mat¬ 
ter halos and galaxies, and it has been remarkably successful 
in describing observations of galaxy clustering (|Seljak| |200 0[ 
|Mo et ar]|2010| for a recent review). In particular, the Halo 
Occupation Distribution (HOD)( [Berlind & Weinberg |2002[ 
[Berlind e t al.||20 03^ and the Conditional Luminosity Function 
(CLF)( Yang et al.|2003^ are the two most widely used models 


of the galaxy-halo connection. These models start from the as¬ 
sumption that halo mass completely determines the galaxy oc¬ 
cupation statistics. In order to populate halos with galaxies, the 
HOD specifies the probability P{N\M) that a halo with mass 
M hosts N galaxies, while the CLF models the mean abun¬ 
dance of galaxies with luminosity L in halos of mass 

M. These two models are interchangeable; integrating the CLF 
over luminosity yields an HOD. Both models have been applied 
extensively to observations in order to study the galaxy-halo 


connection 

^Magliocchetti & Porciani|2003 Zehavi et al. 
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/er, the clustering of halos also exhibits a depen- 
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dence on a 

2005l|Wechsler et al.|2006l|Gao & White|2007l|Li et al.|2008l 


assembly bias. This can be traced back to the fact that halos of 


the same mass in different environments have different assem¬ 
bly histories and cluster differently. Having different assembly 


histories also affects the internal structure of halos([Bullock et al. 


2001l|Wechsler et al.|2002l|Hahn et al.|2007l|Faltenbacher & 

White|2010|l. This, in turn, results in a clustering dependence 

on the stru 
gravitation 
cular veloc 
tion of halo 
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:tural properties of a halo, including the depth of its 
al potential well, characterized by its maximum cir- 
ity Umax. The present work revisits this manifesta- 
assembly bias, and extends it down to smaller scales 
vipc) than has been previously explored. 

emative approach to connecting halos and galaxies 

ce matching (Kravtsov et al.|2004| |Vale & Ostriker 

2004| |Tasitsiomi et al.|2004| Vale & Ostriker|2006| |Conroy & 

Wechsler||2009l IGuo et al.||2010l ISimha et al.||2010l INeistein 

et al.|2010 

[Watson et al.|2012|lRodriguez-Puebla et al.|2012| 

Hearin et a 

.|2012[|Kravtsov|2013[|Saito et al.|2015k In its sim- 

plest form, abundance matching posits a monotonic relationship 
between a property of a galaxy (luminosity, stellar mass) and 


that of a halo (mass, potential well depth). By construction, such 
a relationship preserves the rank ordering of the galaxies and ha¬ 
los. The choice of the observationally-relevant halo property is 
a priori unknown; this uncertainty that can lead to significant 
systematic errors when using halo models to interpret galaxy 
clustering measurements ( [Zentner et al.|201T| . 

A recent parallel research effort has been revisiting the 
standard subdivision of dark matter halos into host/sub-halos, a 
classification that naturally depends on how the halo boundary 
is chosen. While the virial radius is the most commonly chosen 
definition, recent work has demonstrated that the environmen- 
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zel et al.|2014[|Diemer & Kravtsov|20141|Adhikari et al.|2014[ 

■ • high - 

■ ■ - low - Uiiax 

Wetzel & Nagai|2014||More et al.|2015|(. In particular, IWetzell 2.8 

et al. j2014| argue that these environmental effects are due to ._, 



hosts, and therefore get temporarily reclassified as host halos. 
These studies argue that a more physically motivated boundary 
is the “splashback radius” corresponding to the caustic from ma¬ 
terial just reaching its first apocentric passage. One of the chief 
results of the present work is that the halo assembly bias on 
small scales predominantly arises from this mis-classification. 

The outline of this paper is as follows. Sec. 2 summarizes 
the simulations we use in this work. Sec. 3 presents our primary 
results - characterizing the dependence of the clustering of ha¬ 
los on Vjnaxl Sec. 4 explores some of the implications of these 
results. We conclude in Sec. 5. 


2 SIMULATIONS 
We use the Bolshoi (|Klypin et al.|2011a|l and MultiDark sim¬ 


ulations 1 Riebe et al. 


2011 


Prada et al.||20TT] )n in this work; 


the combination of these simulations allows us to span a large 
range in halo mass. We summarize key properties of these sim¬ 
ulations here. Both simulations were run with the Adaptive Re¬ 
finement Tree Code jKravtsov et al.|1997||Gottloeber & Klypin| 
|2008|l assuming a flat ACDM model with density parameters 

= 0.27, riA = 0.73, fit = 0.0469, and ag = 0.82, 
n = 0.95, h = 0.70. The Bolshoi simulation used 2048® parti¬ 
cles in a 250fi“®Mpc box with a force resolution of l/i“®kpc, 
giving a particle mass of 1.35 x while the Mul¬ 

tiDark simulation used 2048® particles in a l/i“®Gpc box 
with a force resolution of 7/i“®kpc giving a particle mass of 
8.721 X lO®/i“®M 0 . Dark matter halos and subhalos are iden¬ 
tified using the ROCKS TAR phase-space, temporal halo finder 
( [Behroozi et al.|[2013d^ and merger trees are constructed us¬ 
ing the CONSISTENT TREES ( [Behroozi et al.|2013c) proce- 
durejf] All of the results we consider here are at 2 : = 0. 

In what follows, we use the virial masses and maximum 
circular velocities (tagged “mvir” and “vmax”) directly from 
the halo catalogues. The halo mass and velocity functions start 
to show incompleteness at lO®®'^/i“®M 0 (~ 200 particles) for 
the Bolshoi simulation, and lO®®/i“®M 0 (~ 100 particles) for 
the MultiDark simulation. Below those masses, the halo mass 
distributions show unphysical drop-offs indicating incomplete 
mass resolution. 

In addition to the standard classification of halos into host 
halos (not within the virial radius of a more massive halo) and 
subhalos (within the virial radius of a more massive halo), we 
further classify host halos into ejected and non-ejected halos. 
Ejected halos, also referred to as “backsplash” halos, are host 
halos whose main progenitor was classified, at some point in 
its merger tree history, as a subhalo. As we discuss below, these 
ejected halos have very different clustering properties compared 
to their non-ejected counterparts. The ejected fraction for the 
Bolshoi simulation is ~ 15.8% at 10^^h~^M q, and drops to 


^ http://www.MultiDark.org 

^ The catalogues used in this paper are publicly available at http; 
//hipacc.ucsc.edu/Bolshoi/MergerTrees.html 
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Figure 1. Distribution of halo mass and maximum circular velocity at 
2 : = 0.0 for halos from the Bolshoi simulation. The blue dots represent 
halos whose observed maximum circular velocity, Vjnax.obs, is greater 
than Umax, while the green dots are the ones with smaller Umax obs 
than Umax ■ The boundary between blue and green dots correspond to 
Umax computed from Eq.[^ 


~ 6.3% at lO®®/i“®M 0 . The lower mass resolution of the Mul¬ 
tiDark simulation prevents us from making this additional sub¬ 
classification. Accordingly, results that rely on this split are re¬ 
stricted to the Bolshoi simulation and mass range. 


3 THE MAXIMUM CIRCULAR VELOCITY 
DEPENDENCE OF HALO CLUSTERING 

In this section we present our primary results. We begin in 
by describing the sample of halos we use throughout the pa¬ 
per, as well as our method for how we categorize halos as hav¬ 
ing above- or below-average circular velocities for their mass. 
In §3.2[ we show the dependence of the clustering of halos on 
Umax. 


3.1 Halo Sample Definitions 

If the internal structure of a dark matter halo of mass Mvir is 
described by an NEW profile ( [Navarro et al.|1997) ) of concen¬ 
tration c, then its maximum circular velocity Vjnax is given by: 


Umax - 0.465MV® jG(g7rAhPcritnm)%®^^^^y~J^YTS) 

( 1 ) 


As shown in ( jKlypin et al.|20lTb| (, the median concentration- 
mass relation c(Mvir) for 2 = 0 Bolshoi halos is well-described 
by: 


log^ijC = —0.0971ogioMvir + 2.148. 


( 2 ) 


Eor every halo in the Bolshoi and MultiDark catalogs, we use 
its tabulated Mvir to compute c(Mvir), and then use the values 
c together with Eq.[^to compute Umax for every halo. We will 
henceforth refer to halos with Umax < Umax (Mvir) as “low- 
Umax halos”, and halos with Umax > Lmax(Mvir) as “high- 
Umax halos”. Thus a halo’s high- or low-Umax designation refers 
to whether its true Umax value in the simulation is above- or 
below-average for its mass. 
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Fig. [^illustrates the distribution of Bolshoi halos as a func¬ 
tion of Mvir and l/max- High-V^ax halos are shown in blue, 
low-Vinax halos in green. The dividing line between the two 
samples is defined by Eq.[T] We find that our analytical approx¬ 
imation for l^max(Afvir) gives a good description of the true 
median: for any fixed value of Mvir, the high-Finax and low- 
Fmax subsamples have very similar numbers of objects. 

Note that Fig. [T] shows a sharp lower bound on the value 
of Vmax at a given Mvir, but no sharp upper bound. This is ul¬ 
timately due to the halo mass definition. The circular velocity 
14ir at the virial radius of any halo is 14ir(F?vir) = Fvir = 
f GMyir/Rvir- Since the value of Knax tabulated in the halo 
catalog is computed as the maximum value of Mir over the en¬ 
tire profile of the halo, formally Vmax cannot exceed Mir. This 
manifests as the sharp lower bound seen in Fig.[T] 


3.2 Halo Bias 

In this section we present our primary results for the cluster¬ 
ing properties of halos as a function of Mvir and Vjnax. Clus¬ 
tering strength is quantified by the two-point correlation func¬ 
tion, ^(r). In all that follows, we will use to denote the 

auto-correlation of the dark matter density field with itself, and 
Chm(f) to denote the cross-correlation between a sample of ha¬ 
los and the underlying density field. 

Halos are biased tracers of the dark matter density field. We 
denote this bias as which is in general a function of spatial 
separation. We define halo bias as 

bh(r) = Chm(r)/(^mm(?'). 

On sufficiently large scales halo bias is approximately linear, 
and bh(r) approaches a constant value 

In order to measure the bias of a sample of simulated ha¬ 
los for both Bolshoi and MultiDark simulations, we estimate 
Cmm(r) and 5 hm(f) using a random down-sampling of 10 ® dark 
matter particles. For a given sample of halos, we estimate the 
value of fejj" exhibited by the sample as follows: 

&lin = ^ ^ (Chm(fi)/Cmni(o))- (3) 

i 

In Eq.j^ the sum is performed over A^bin = 20 separation bins 
n linearly spaced from 10/i~^Mpc to 20h“^Mpc. 

In order to study the mass-dependence of halo bias, we bin 
our halos into a sequence of Mvir bins chosen such that there 
are the same numbers of halos in each bin. For MultiDark, we 
select 2 x 10® halos for each bin; for Bolshoi we use 25000 
halos per bin. The halos in each mass bin are categorized as 
high-Mnax or low-Miax according to the method described in 
We start using halos from the MultiDark simulation for 
Mvir > lO^^-^/i-^M 0 . 

The top panel of Fig. [^ shows as a function of Mvir; 
results for high-Mnax halos are shown in blue, low-Mnax ha¬ 
los in green. At the low-mass end, linear bias is a weak func¬ 
tion of Mvir; for Mvir > Mcoii « IO^^'^Mq/Zi, where Mcoii 
is a characteristic mass scale for clustering corresponding to 
O'(Mcoii, 2 ) = Sc ~ 1.69, we see that increases sharply 
with Mvir. Thus the basic shape of each curve in Fig. [^ is in 
accord with theoretical expectations from the peak-background 
split ( |Sheth & Tormen||1999|> fo rmalism and Press-Schechter 
theory ( [Press & Schechter||1974| l. By comparing the blue and 
green curves in Fig.j^we can see that linear bias has significant 
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Figure 2. Upper panel: Linear bias at z = 0 as a function of halo 
mass from the Bolshoi and the MultiDark simulations. The vertical 
dashed line indicates the upper/lower limit of halo mass for the Bol- 
shoi/MultiDark simulations respectively. The blue lines correspond to 
a linear bias for high-Vinax halos, while the green lines correspond 
to low-Mnax halos. The solid lines correspond to the full halo sample 
while the dashed line corresponds to the sample where ejected halos ai'e 
removed. Lower panel: Ratio of linear biases between high-Vjnax and 
low-Vjnax samples from the Bolshoi and the MultiDark simulations. 
The clustering difference increases with decreasing Mvir, and exceeds 

30% for halos of Milky Way mass Mvir ~ fO^^MQ/h. 


dependence upon Knax for halos of the same mass. This depen¬ 
dence is most readily seen in the bottom panel, which shows the 
ratio of of high-Umax samples divided by fe®" of low-Mnax 
samples. Thus in the bottom panel of Fig.j^ vertical axis values 
exceeding unity correspond to masses where high-Mnax halos 
cluster more strongly relative to their low-Vjnax counterparts. 

At the low-mass end, high-Mnax halos are more strongly 
clustered than low-Vjnax halos of the same mass. The clus¬ 
tering difference increases with decreasing Mvir, and exceeds 
30% for halos of Milky Way mass Mvir ~ IO^^Mq/Zi. At 
the high-mass end, the trend reverses, and the overall magni¬ 
tude is weaker. These results are consistent with lWechsler et ah] 
j2006| l, who find that the same trends hold when the secondary 
halo property is NFW concentration, rather than Vjnax. This 
agreement is to be expected: insofar as the halo profile is well- 
approximated by an NFW profile, Vjnax is entirely determined 
by concentration (see Eq. [TJ. Note that there is a drop in the 
ratio of the linear biases on the low-mass end. As we see no 
physical reason for this drop, we consider this to be a resolu¬ 
tion effect and suggest that the completeness requirements for 
two-halo-property-dependent clustering are significantly more 
stringent relative to the requirements demanded by the need for 
a complete halo mass/velocity function. 

Next, we study the scale-dependence of halo bias on small 
scales for high-Vjnax and low-Vjnax halos. We wish to parse 
the novel, small-scale effects from the well-known large-scale 
effects. So for each sample of high- and low-Vjnax halos, we 
compute the following quantity: 
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&h(r|Mvir; Knax) = 6h(r|Mvir; Knax)/&lin (Mvir; Kiax)- 

(4) 

Thus for any sample of halos, as r > lOMpc/li, we have 
bh(r) —> 1, by construction. 

The bottom panel of Fig. shows the ratio of bh{r) of 
high-Knax samples divided by hh{r) of low-Vmax samples for 
several mass bins. The first three mass bins labeled in the figure, 
Mvir = 10^^'^42 0 ,i 2 . 2 ^-i]y|-^^ from the Bolshoi simula¬ 
tion, and the last two mass bins, Mvir = 

from the MultiDark simulation. High-Vmax halos cluster more 
strongly compared to their low-Knax counterparts at l/i“^Mpc. 
This scale-dependent feature becomes stronger with decreas¬ 
ing Mvir and exceeds 40% for halos of Milky Way mass and 
reaches 60% at Mvir ~ 

Up until present, we have used both host halos and ejected 
halos to compute halo biases. Both types of halos are identified 
as distinct halos at 2 = 0. Ejected halos, however, are halos 
which were identified as part of more massive halos at one or 
more occasions in the past, but were ejected and now exist as 
a host halo at 2 = 0. Those ejected halos tend to exist around 
more massive halos jWetzel et al.|2014^|Wang et al.|2009] e.g.). 
Therefore, the effect on scale-dependent biases may be caused 
by those ejected halos. 

To test this ejected halo hypothesis, we compute halo- 
matter cross correlation functions after first excluding the sub¬ 
population of ejected halos. Our results for the linear regime are 
shown as the dashed curves in Fig.[^ The relative difference in 
the linear bias between high-Knax and low-Knax halos is sup¬ 
pressed to 25% for halos of Milky Way mass. This suppression 
due to excluding the ejected halos is consistent with the results 
presented in | Wang et al.| ( |200^ . 

In the bottom panel of Fig. we show the scale- 
dependence of assembly bias for non-ejected halosj^Once the 
ejected halos have been removed, the scale-dependent feature 
of assembly bias is greatly reduced. This implies an intimate 
connection between the scale-dependence of assembly bias and 
subhalo back-splashing (see ^for further discussion). 


4 OBSERVATIONAL CONSEQUENCES 


We now consider possible observational consequences of the re¬ 
sults of the previous section. In order to do this, the key first step 
is to relate an observable property of a galaxy (luminosity, stel¬ 
lar mass) to an intrinsic property of its host halo (mass, circular 
velocity). As one might infer from above (and we demonstrate 
below), different choices for the latter can result in significant 
differences for different observables. 

In order to be explicit, we use the abundance matching 


Tasit- 


siomi et al.||2004[ |Vale & Ostriker||2006 

Conroy & Wechsler 

20091 |Guo et al.|20101|Simha et al.bolO 

Neistein et al.|2010| 

Watson et al.||2012[ |Rodrlguez-Puebla et al.||2012t |Kravtsov 


the mass or circular velocity of host halos. We implement this 
by splitting the halo catalog into a series of bins with constant 


® We remind the reader that due to the mass resolution of the MultiDark 
simulation, fig.j^only shows results for Bolshoi halos. 



Figure 3. Top panel: Ratio of halo-matter cross correlation functions 
between high/low-VKax halos from the Bolshoi simulation and the 
MultiDark simulation at 2 = 0, normalized by their linear biases. Each 
line corresponds to different halo mass bins labeled in the plots. The 
top three lines that correspond to low mass halos are computed from 
halos in the Bolshoi simulation, while the bottom lines are from the 
MultiDark simulation. Those plots show that high-VKax halos clus¬ 
ter more strongly than low-Knax halos at lh“^Mpc and the relative 
scale-dependence between those subsamples increases smoothly with 
decreasing halo mass. Bottom panel: The same figure as the top panel 
without ejected halos only from the Bolshoi simulation. As can be 
seen by comparing these results to those in the top panel, the Knax- 
dependence of halo bias on small scales is dramatically reduced by ex¬ 
cluding ejected subhalos. This implies an intimate connection between 
the scale-dependence of assembly bias and subhalo back-splashing. 


number density (=1.6 x 10“^(/i“^Mpc)“®), rank ordering ei¬ 
ther by mass or circular velocity. We label each bin by its corre¬ 
sponding stellar mass, computed from the stellar-to-halo mass 
relation of ( [Behroozi et al.|2013at . 

Note that when we rank order based on circular veloc¬ 
ity, there is the possibility that the mean halo masses of these 
bins could differ from what we obtain after rank ordering by 
halo mass. We explicitly check this and find that the mean halo 
masses for both cases agree to ~ 99.6%, allowing us to consis¬ 
tently compare samples of mock central galaxies with the same 
stellar mass, but where the stellar mass is statistically regulated 
hy either Mvir or Knax. 

We find only a relatively minor difference in the large-scale 
clustering of the two samples of mock central galaxies. At fixed 
stellar mass, the linear bias of samples selected by their circular 
velocity are ~ 5% higher than samples selected by halo mass. 
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Figure 4. Ratio of halo-matter cross correlation functions between 
Mvir-based and ^max-based abundance matching samples with ejected 
halos (top) and without ejected halos (bottom). The ratios ai‘e nor¬ 
malized by their linear biases. Each line corresponds to different stel¬ 
lar mass bins labeled in the plots. With ejected halos, the halo biases 
computed from the Vmax-selected samples have very different scale- 
dependence than the ones from the Myir-selected samples. By remov¬ 
ing those ejected halos, the difference is reduced by threefold. 


This decreases to ~ 2% if we remove ejected halos from both 
samples. 

We study the scale-dependence of the clustering of our 
mock central galaxies in Fig. which is directly analogous to 
Fig. a only here we have use the abundance matching tech¬ 
nique described above to illustrate how our “halo-level” results 
may manifest in observed galaxy populations. Again we see a 
clear scale-dependence of the clustering signal, with a maxi¬ 
mum difference of ~ 15% at ~ lh“'^Mpc. These differences 
go down to ~ 5% after removing ejected subhalos, again re¬ 
flecting the intimate connection between scale-dependent as¬ 
sembly bias and subhalo back-splashing. 


5 DISCUSSION 

For halos Mvir ^ Mcoii ~ we have shown 

that the linear bias of low-Vjnax halos is larger relative to high- 
Vniax halos of the same mass. As shown in |Dalal et al.| ( |2008t , 
this phenomenon is nicely explained in terms of the statistics 
of fluctuations in a Gaussian random field. Consider two halos 
with the same present-day mass, but with different concentra¬ 
tion. Both halos originate from a fluctuation of the same peak 
height, but with different peak curvature: the high-concentration 


(high-Umax) halo has a sharper peak than the low-concentration 
(low-Vniax) halo. |Dalal et al.|(2008| l showed that a generic pre¬ 
diction of Extended Press Schechter theory (EPS) with a con¬ 
figuration space filter is that low-curvature peaks cluster more 
strongly relative to high-curvature peaks of the same height. In 
closely related work, |Zentner| ( [2007) l used EPS with a conhgu- 
ration space filter to show that for a pair of halos of the same 
peak height, the early-forming halo should reside in a denser 
large-scale environment than the late-forming one[^ 

A critical assumption underlying these EPS predictions is 
that a halo is the dominant peak in its large-scale environment. 
This is a well-founded assumption at the high-mass end, and we 
see that the predictions are in good agreement with simulations 
in the Mvir > Mcoii regime jPalal et al.|2008^ . The situation is 
quite different when Mvir < Mcoii. We have confirmed previ¬ 
ous results (e.g., |Wechsler et al.|2006^ that large-scale assembly 
bias changes sign and strengthens for lower-mass halos. This 
is in stark contrast to the EPS model described above, which 
makes the same prediction for assembly bias regardless of halo 
mass. 

Thus EPS succeeds and fails in precisely the regimes 
where we expect. Lower-mass halos are strongly influenced by 
the tidal field in which they evolve iHahn et al.||2007 Wang 


let al.|[20TT] |Shi et alT||2015| |Hahn et al.||2009| [Hearin et al. 


|2015^ ; the EPS assumption that the halo dominates its environ¬ 
ment breaks down catastrophically, and in this regime nonlinear 
evolution governs assembly bias. On the other hand, high-mass 
halos do dominate their tidal environment; the EPS assumption 
holds good, and we can understand assembly bias as naturally 
arising from the statistics of Gaussian fluctuations. 

Our results on the scale-dependence of assembly bias are 
also consistent with this picture. Eirst, we remind the reader that 
halo bias for high-Vjnax halos shows non-trivial scale depen¬ 
dence with a pronounced bump at ~ 1 — 2Mpc//i compared 
to the bias for low-Umax halos. This scale-dependent feature for 
high-Unax halos becomes 60% larger compared to low-Vmax 
halos at Mvir ~ IO^^'^Mq/Zi. This feature, however, is re¬ 
moved by excluding the ejected halos, implying that this special 
sub-population is responsible for the scale-dependent bump. 

This scale-dependence has a simple interpretation in terms 
of subhalo back-splashing. Eirst, ejected halos are physically 
associated with the more massive halo from which they were 
ejected. The clustering of ejected halos is therefore largely de¬ 
termined by this associated massive halo, much like the clus¬ 
tering of present-day subhalos is determined by their host halo. 
Second, as ejected halos pass near and inside a massive halo, 
their physical growth is arrested, and many such halos even ex¬ 
perience substantial mass loss ( |Wang et al.|2009||Wetzel et al.| 
|2014| ). This arrested development has a greater impact on the 
outer layers of the halo, so that the ejected halo’s mass is signif¬ 
icantly more affected than its circular velocity jBehroozi et al.| 
|2013b^ . Putting these two effects together, we should naturally 
expect the outskirts of massive groups and clusters to be prefer¬ 
entially populated with low-mass halos that have above-average 
values of Umax for their mass. This manifests in the scale- 
dependent bump shown in Figures[^&|^ 

Recent advances in our understanding of halo growth 
sheds further light on the above results. As shown in[D 


|& Kravts^ p014| l; |Adhikari et al.| ( |2014^ ; | Wetzel & Nagai 


^ See Chapter IX, Section D. 
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( |2014^ ; [More et al.| ( |2015[ l, the natural physical boundary of 
a dark matter is the so-called “splashback radius”, which is 
the radius where accreted matter reaches its first apocenter af¬ 
ter turnaround, and is roughly 2 — 3-Rvir. As shown in |Wet-| 
|zel et ah] \1Q\A) , the halos of massive groups and clusters 
(Afvir ^ 'W^Mq / h, Rvir > 500kpc//i) are surrounded by 
a large fraction of ejected halos. For ejected halos with Mvir ~ 
halos, the host-centric spatial distribution of the 
ejected population peaks at r ~ 1.25i?200m ~ 1.5Mpc/fe (see 
Figs. 2 and 3 in [Wetzel et al.| ( |2014^ ). The bump feature we 
find at 1 — 2Mpc//i is therefore in quantitative agreement with 
the [Wetzel et al.|j2014^ ; |Wang et al^ ( |2007| ) results: this bump 
occurs at the same physical scale that we would expect if the 
clustering of the ejected population is largely determined by the 
host to which the halos are ultimately bound. Note that the as¬ 
sembly bias on large scales, r > 10Mpc//i, is not dominated by 
the ejected halos, which is consistent with [Wang et al.| ( |2007[ >. 
Without the ejected halos, the relative difference in the linear 
bias between high-Vmax and low-Vinax halos remains 25% for 
halos of Milky Way mass. 

In order to explore possible observational consequences 
of our findings, we use abundance matching relating a stellar 
mass of central galaxies to either mass or circular velocity of 
host halos. Using different intrinsic properties for host halos re¬ 
sults in significant differences in clustering signals of central 
galaxies. On large scales, the linear bias of samples selected 
by their circular velocities are ~ 5% higher than samples se¬ 
lected by halo masses at fixed stellar mass. Without the ejected 
halos, this differences is reduced to ~ 2%. On small scales 
(r < 10/i“^Mpc), samples selected by their circular velocities 
exhibit the scale dependence with a bump at 1 — 2 /i“'^Mpc com¬ 
pared to samples selected by halo masses. This scale-dependent 
bump for Vmax—selected samples becomes ~ 15% with the 
ejected halos and ~ 5% without the ejected halos. As these 
effects are roughly as large as existing SDSS clustering mea¬ 
surements on these scales, this raises the possibility that clus¬ 
tering measurements can be used to determine which host halo 
property is the true statistical regulator of the stellar mass of the 
central galaxy residing in the halo. 


6 SUMMARY 

We conclude the paper with an overview of our primary results: 

(i) At fixed mass Mvir, the large-scale bias of halos exhibits 
significant residual dependence on potential well depth Vinax. 
At the low-mass end, high-Umax halos cluster more strongly 
than their low-Umax counterparts. At the high-mass end, this 
trend reverses, and is generally weaker, with the transition oc¬ 
curring at Mcoii ~ lO^^'®M 0 //i. Our results are quantitatively 
consistent with previous studies of large-scale halo assembly 
bias. 

(ii) We show that assembly bias exhibits complex scale- 
dependence. The Umax—dependence of halo clustering shows 
a pronounced “bump” on scales 500kpc//i ^ r < 5Mpc//i. 
This scale-dependence is itself mass-dependent: the bump fea¬ 
ture is strongest for low-mass halos and vanishes for halos with 
Mvir > Mcoll. 

(hi) The scale-dependence of assembly bias can primarily be 
attributed to a special sub-population of ejected subhalos, which 
experience arrested mass-growth before and after being ejected 


from a higher-mass host. If this special population is excluded 
from the halo sample, the strength of small-scale assembly bias 
is limited to < 5% for all masses Mvu Mq/Ii. 
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